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^! Abstract 

D 

Long-distance contributions to the D^-D^ mixing parameters x and y are evaluated using lat- 
C^ ■ est data on hadronic D^ decays. In particular, we take on two-body D — )■ PP and VP decays 

\^ • to evaluate the contributions of two-body intermediate states because they account for ~ 50% of 

^P ' hadronic D^ decays. Use of the diagrammatic approach has been made to estimate yet-observed 

decay modes. We find that y is of order a few xlO^^ and x of order 10~^ from hadronic PP and 
lO ' VP modes. These are in good agreement with the latest direct measurement of D^-D^ mixing 

^— s ' parameters using the D^ — )• i^^vr+vr" and KsK^K~ decays by BaBar. We estimate the contri- 

bution to y from the VV modes using the factorization model and comment on the single-particle 
resonance effects and contributions from other two-body modes involving even-parity states. 



I. INTRODUCTION 

It is well-known that short-distance contributions to the D^-D^ mixing parameters x and y to 
be defined below are very small, of order 10~^, jawing to the Glashow-Iliopoulos-Maiani (GIM) 
suppression and the doubly Cabibbo suppression ^i, i^].^ Since the mixing effects and CP violation 
in the neutral charmed meson system are expected to be very small compared to kaon and B 
mesons, it is difficult to observe them experimentally. Nevertheless, BaBar 5(], Belle [6| and CDF 
[7|] have provided compelling evidence for D^-D^ mixing in the past few years. The current world 
averages for the CP allowed case are [8(] 

X = (0.98l°J^)% , y = (0.83 ± 0.16)% . (1) 

However, the analyses that have reported evidence for mixing have not been able to provide direct 
measurements of x and y. A most recent BaBar experiment using the D^ — )• Kgn^n^ and D^ — >■ 
KgK^K^ decays measured the mixing parameters x and y directly, with the results [9[ 

x = (1.6 ±2.3 ±1.2 ±0.8) X 10"^ , y = (5.7 ± 2.0 ± 1.3 ± 0.7) x 10"^ , (2) 

where the last error comes from the amplitude models used in the analysis. These results are 
consistent with the previous Belle measurements from the i^^vr^vr^ mode alone [10|: 

X = (8.0 ± 2.9l:Ut\i) X 10-=^ , y = (3.3 ± 2.4lO|t[J;^) x 10"=^ . (3) 

Therefore, the new BaBar measurement favors lower values for x than for y and lower x and y 
values than the aforementioned world averages. At any rate, the observed D^-D^ mixing is several 
orders of magnitude larger than what is expected from the short-distance contributions, and is 
evidently dominated by long-distance processes. 

Since the long-distance effects are non-perturbative in nature, it is conceivable that it will be 
very difficult to have a reliable estimate of the charm mixing parameters, especially in view of the 
fact that the charm quark is not heavy enough for a sensible heavy quark expansion and not light 
enough for the application of chiral perturbation theory. In general, the long-distance contributions 
are estimated in either of the two approaches: inclusive and exclusive. The "inclusive" approach 
relies on the heavy quark expansion dictated by the parameter l/rric (for a recent study, see 



111]). In the so-called "exclusive" approach, on the other hand, one sums over contributions from 
intermediate hadronic states. Since there are cancellations among states within a given flavor SU(3) 
multiplet, as first noticed in J12l. Il3l|. one needs to know the contribution of each state with high 
precision in order to have a trustworthy estimate |lj| . It can be shown that the mixing parameters 
X and y vanish in the SU(3) limit. In the exclusive scenario, this means that the cancellation of the 
Cabibbo-favored (CF) and doubly-Cabibbo-suppressed (DCS) decays with the contributions from 
singly-Cabibbo-suppressed (SCS) transitions is perfect in the limit of SU(3) symmetry. In other 
words, nonzero values of x and y comes from the breakdown of flavor SU(3) symmetry. In the 
absence of sufficiently precise data on many decays rates and on strong phases, the authors of [ij] 



^ For a tabulation of predictions on the mixing parameters x and y within and beyond the standard model, 
see 0,3- 



computed the long-distance contributions to AF by considering SU(3) breaking from final-state 
phase space differences and neglecting SU(3) violation in the decay amplitudes. They found that 
the phase space effects alone could produce sufficient SU(3) breaking to induce y ~ 10"^ and that 
large effects in y appeared in decays to final states close to the D threshold. 

In the past few years, rich data on hadronic D — t- PP, VP decays have been accumulated from 
various experiments with improved accuracy. (Throughout the paper, we use P, V, S, A, and T 
to denote pseudoscalar, vector, scalar, axial-vector, and tensor mesons, respectively.) Especially, 
there are new CLEO measurements for the PP modes with better precision, many of them having 



experimental errors less than the present world averages 15(]. These data allow us to make a 
sensible determination of the mixing parameters x and y from the PP and VP channels without 
relying on model assumptions. ^ There are some channels that have not been measured: six SCS 
VP modes and many of DCS PP and VP decays. We will employ the diagrammatic approach to 
give inputs for those unmeasured channels. 

Recently, we have studied the two-body hadronic charmed meson decays, includin g al l the PP, 



VP, SP, AP and TP modes, within the diagrammatic and factorization approaches [19|, |20||. The 
best-fitted values extracted from the CF decay modes in the diagrammatic approach have been 
used to predict the branching fractions of the SCS and DCS modes for the D — t- PP and D — t- VP 
decays. This approach enables us to estimate the mixing parameters x and y. 

The layout of this work is as follows. In Section II, we first review the method of computing 
y from hadronic decay branching fractions, and then consider various two-body intermediate-state 
contributions to y. In Section III, we review the dispersion relation for x, and compute PP and VP 
mode contributions to x, followed by a brief discussion on single particle effects. We summarize 
our findings in Section IV. Our Appendix A compares different state normalization conventions 
commonly seen in the literature. Appendix B discusses the relative strong phase between the 
D^ — 7- K^TT^ and K^tt^ modes. 

II. DECAY WIDTH DIFFERENCE 

The neutral D meson mass eigenstates are defined in terms of flavor eigenstates as 

\Di,2)=p\D'')±q\D'). (4) 

Using the mass and width matrices, the ratio q/p reads 



where the convention CP\D^) = \D ) has been made. The parameters x and y are defined as 

Am mi - 771,2 AF Fi - F2 ,„. 

x^ — = ^^—, y^^ = ^^- (6) 



For early attempts to estimate the long-distance contributions to a; or y from two-body states, see 



16-1 



Since CP violation in both D mixing and decays is expected to be small within the standard 
model and in most new physics scenarios, we therefore define the CP eigenstates 

|Z)±) = ^(l^°> ± \d')) , with CP\D±) = ±\D±) . (7) 

Hence Di ~ D^ and D2 ~ D^. In perturbation theory, the off-diagonal mass and width matrix 
elements are given by [21i]'^ 

M-M =^(I)»|fl„,|lf)+^v'<"°l^-l")<"l"-l^). (8) 

2/12 2mD 2mD'^M mo-En + ie ' ^ '^ 

where AA is an appropriate normalization factor for the intermediate state \n); for example, A^ = 
2En for a one-particle intermediate state. Using the relation 

V — - indiniD - En) , (9) 



niD - En + ie mo - En 

with V denoting the principal value prescription, we have 



Ar = ^ Y. Jf [{D^\Hn,\n){n\H^\^) + 0'\H^\n){n\H^\D'^)\ (27r)5(mD - En) . (10) 

Therefore, Am, and AF are induced by off-shell and on-shell intermediate states, respectively. 
The parameter y has the expression 



2r 
2r 



^^p, {\{D+\H^\n)\-' - \{D^\H^\n)\^) , (11) 



where pn is a phase-space factor. For example, pn = Pc/(8vrm^) for the D — )• PP decays, with pc 
being the center-of-mass momentum of either meson in the final state. Defining CP\f) = 7ycp|/)) 
y can be recast to [ij] 

y = :^Y.Pn^Cp{n){{D'^\H^\n){n\H.^\D^) + {D^\H^\n){n\H^\D'^)) 



= E ^CKM(n)r/cp(n) cos 5n^B{D^ -^ n)B{D^ ^ n) , (12) 

n 

where 5n is the strong phase difference between the D^ — )• n and D^ — )• n amplitudes and r/cKM = 
(—1)""* with n^ being the number of s and s quarks in the final state. The factor rycp = ±1 is 
well-defined because |/) and |/) are in the same SU{3) multiplet. This factor is the same for the 
entire multiplet. 



^ The expressions of the D'^-D matrix element given in the literature are often very confusing as they are 



not dimensionally consistent. This issue is discussed in Appendix A. In practice, we do not use Eq. (|T0] 
to compute the long-distance contributions to Atti and AF. Rather, we employ Eq. ([12]) and Eq. (l36l 
[or Eq. p7p ] to evaluate the parameters y and x, respectively. They are free of any ambiguities. 



A. PP 

Since CP\Ti^) = — Ivr'^) and likewise for r/, ry', we will choose the convention such that CP\K^) = 
— \K~) and CP\K^) = —1-^"). Because under the CP transformation 

CP\MiM2) = ?7cp(Mi)77cp(M2)(-l)^|MiM2) = r/cp(MiM2)|MiM2) , (13) 

it is clear that r]cp{PP) = 1 for decays into two pseudoscalar mesons. The parameter y arising 
from the PP states is 

ypp = BiTT+TT-) + fi(7rV°) + e(7r%) + B{7r%') + B{r]r]) + ^(W) + B{K+K~) + B{K^K°) 

-2 cos 6 KO^^B{K^r])B{K^r]) - 2 cos 6 ko r,' ^ B{KS)B{KS) ■ (14) 

To see that ypp vanishes in the SU(3) limit, we should work on the SU(3) singlet state ryo and 
octet states 7T,K,rjs. The octet states have the same masses when SU(3) symmetry is exact. We 



2 cos 6K+n~ J B{K-TT+)B{K+TT-) - 2cos6ko^oJ B{KOttO)B{KOtt 



first write down the general quark-graph amplitudes (see J19l ] for details): 



A{D^ ^ K-7r+) = V:,VUT + E) , A(Z?o ^ i^V) = i=T/* K^(C - E) , 

A{D^ ^ K%8) = -^V:,VUC - E) , AiD"" ^ :^%) = -^V:,Vud{C + 2E) , 






^(Z)0 ^ ^+^-) ^ KdKd(T' + ii;') , ^(Z)° ^ 7r°7r°) = i=y;^Krf(C' - i?') , 



^(1)0 ^ TrOryo) = -^KdKdii^' + ;7|^-KsC" , 

A{D^ ^ 7,8%) = ^l/;,Kd(C' + E') + ^y;^K.(-2C' + AE') , 

A{D^ ^ rygryo) = -^V:,VUC' + E') + -i=y;,Ks(C" + 4i?') , 

^(1)0 ^ ^+^-) ^ k:ks(t' + E') , A{D^ ^ ^0;^°) = v:^VudK + K:Ksi?d , 

A(D0 ^ K+vr-) = y;,Ks(T" + E") , A{D^ ^ K^n') = ^V^^VusiC" - E") , 

A{D' ^ i^^s) = -^V:,Vus{C" - E") , A{D' ^ i^^o) = ^y;,Ks(C" + 2E") , (15) 

where T, C, E are color- allowed, color-suppressed and W^-exchange amplitudes, respectively. We 
have followed the conventional practice to denote the primed amplitudes for the SCS modes and 
double-primed amplitudes for the DCS decays. In the SU(3) limit, the primed and unprimed 
amplitudes should be the same, the strong phase 5„ vanishes, and the D^ — )• K^K decay is 
prohibited. It is easily seen that perfect cancellation occurs among the SU(3) octet final states, 
as well as among the decay modes T^^rjo, r}^r]Q, K ryo and K'^ryo involving the SU(3) singlet ryo. 
Therefore, ypp indeed vanishes in the SU(3) limit. 

Flavor SU(3) symmetry breaking occurs in both the decay matrix elements and in the final-state 
phase space. Phase space is the only source of SU(3) violation considered in the previous analysis 



of [ij], given the paucity of data at that time. Since ah the data of D ^ PP are now available 
except for three of the DCS modes, we can have a much more accurate estimate of ypp directly 
from the data. For the yet to be measured DCS modes, we can rely on the relations, for example, 
T{D^ — )• K^7r^)/T{D^ — )• K tt^) = tan^ 9c that has been tested in the measurement of the quantity 

P(n') = r(ZJ° ^ Ksn^) - T{D' ^ i^z^vrO) 
The prediction R{D^) = 2tan^ Oq = 0.107 agrees quite well with the experimental value of 0.108 it 



0.025 ± 0.024 by CLEO [23]. 



Prom the model-independent analysis in the diagrammatic approach, it has been observed that 
sizable violation of flavor SU(3) symmetry occurs in some of SCS modes. The most noticeable 
example is the ratio R = r{D^ — t- K^K^)/T{D^ — t- vr^vr^) ~ 2.8. If the SU(3) symmetry breaking 
manifested only in the phase space, one would obtain R = 0.86 and hence the K^K^ production 
rate should be smaller than the vr^vr^ one. This is in sharp disagreement with experiment. We 
have shown in 19|] that in addition to the SU(3) breaking effect in the spectator amplitudes, the 
long-distance resonant contribution through the nearby resonance /o(1710) can naturally explain 
why D^ decays more copiously to K^K~ than vr+vr" through the VK-exchange topology. This has 
to do with the dominance of the scalar glueball content of /o(1710) and the chiral-suppression effect 
in the decay of a scalar glueball into two pseudoscalar mesons. The same final-state interaction 
(FSI) also explains the measured rate of D^ — )• K^K^ even though its amplitude vanishes in the 
SU(3) hmit. 

CLEO has measured the relative strong phase between D^ — t- K^tt^ and D'^ — t- K^ti^ to be 
cos 5 = l.OSl^Qj^y lb 0.06 [24]. (See Appendix B for an estimate of this strong phase.) It is thus 
plausible to assume cos 5^ = 1 for all n = PP. From Eq. ()14p and the data in Table HI we obtain 

ypp = (1.128 ± 0.038)% - (1.042 ± 0.017)% = (0.086 ± 0.041)% . (17) 

B. VP 

The neutral vector mesons p^,uj,(l) are CP eigenstates with CP = +. It is thus convenient to 
define CP\V) = \V) for the vector mesons in the same SU{3) multiplet. It follows from Eq. (fT3]) that 
r]cp(VP) = +1 for decays into one vector meson and one pseudoscalar meson. There are more decay 
modes available for the VP final states, namely, V1P2 and P1V2. There are totally 30 VP channels 
(8 for CF, 14 for SCS and 8 for DCS), to be compared with 16 PP channels. Because the decay 
constant of the vector meson /y, typically of order 210 MeV, is much larger than /p, many VP 
modes have rates greater than the PP ones. For example, B{K~ p~^) ~ 11% >> B{K~tt~^) ~ 4% and 
B{7r~ p'^) ~ 1% » B{7r~^7r~) ~ 1.5 x 10"^. It is thus anticipated that the VP mode contributions 
to y will dominate over ypp. 

Note that the decay amplitude of the DCS mode is not simply related to the corresponding CF 
one by replacing the CKM matrix elements V*gVud with V*^Vus- For example, 

A{dO -^ K*+7r-) _ V;^Vus Tp + Ev 



A(DO ^ K*-7T+) V*,Vud Tv + Ep' 



(18) 



TABLE I: Experimental branching fractions for Cabibbo- favored (in units of %), singly-Cabibbo- 
suppressed (in units of 10~^) and doubly-Cabibbo-suppressed decays (in units of 10"'^) of D^ — )■ 
PP, VP. Data are taken from 15| for D — )• PP and from [23] for D — ?• VP. The channels with the 
superscript * have not been measured yet. For them, we use the fitted branching fractions obtained 
from the diagrammatic approach [19]. For D — )• VP decays, the fitted branching fractions are those 
obtained from solution (A,Al) (upper row) and (S,S1) (lower row) in [l9[]. 



Mode 



B, 



cxp 



Mode 



B. 



exp 



Mode 



B, 



cxp 



K rj 



^OyrO 



7rO?7' 



rjT] 
rjrj' 
K+K- 



K+ 



vr 






(3.91 ± 0.08)% 
(2.38 ± 0.09)% 



(0.96±0.UDj7o 

(1.90 ±0.11)% 

1.45 ±0.05) X 10-3 

0.81 ±0.05) X 10-3 
0.68 ± 0.07) X 10-3 
0.91 ±0.13) X 10-3 

1.67 ±0.18) X 10-3 

1.05 ±0.26) X 10-3 

4.07 ±0.10) X 10-3 

0.64 ±0.08) X 10-3 

1.48 ± 0.07) X lO-'' 



0.67 ±0.02) X 10- 



0.28 ±0.02) X 10 



0.55 ± 0.03) X 10" 



-4* 



-4* 



'4* 



K-p+ 

K T] 

7r-p+ 



7rV° 



K+K*- 
K-K*+ 

Km*"" 



tt'^oj 



(5.91 ±0.39)% 
(2.82 ±0.35)% 

(10.8 ±0.7)% 

(1.54 ±0.12)% 

(0.96 ± 0.30)% 

(0.012 ± 0.003)% 

(2.26 ± 0.40)% 
(0.868 ±0. 



(4.97 ±0.23) X 10-3 

(9.8 ±0.4) X 10-3 
(3.73 ±0.22) X 10-3 
(1.53 ±0.15) X 10-3 

(4.41 ±0.21) X 10-3 

(0.05 ±0.06) X 10-3* 

(0.29 ± 0.22) X 10-3* 

(0.29 ± 0.22) X 10-3* 

(0.05 ± 0.06) X 10-3* 

(0.10 ±0.18) X 10-3* 

(1.01 ±0.18) X 10-3* 



TjU) 



rj'ui 



r]p'' 



7]'p^ 



K*+Tr- 



P 



K+ 



p^K^ 
6K^ 



(1.24 ±0.12) X 10-3 
(2.21 ±0.23) X 10-3 
' (0.07 ± 0.02) X 10-3* 
(0.15 ±0.01) X 10-3* 
(0.14 ±0.05) X 10-3 
(1.11 ±0.86) X 10-3* 
(1.17 ±0.34) X 10-3* 
(0.14 ±0.02) X 10-3* 
(0.26 ±0.02) X 10-3* 



(3.54l{; 



80 \ 
05^ 



X 10' 



(0.54 ±0.18) X 10-^* 
(0.74 ±0.17) X 10-^* 
(1.45 ±0.17) X 10-^* 
(1.91 ±0.21) X 10-^* 
(0.91 ±0.51) X 10-^* 
(0.63 ±0.19) X 10-^* 
(0.33 ± 0.08) X 10-^* 
(0.28 ± 0.05) X 10-^* 
(0.0040 ± 0.0006) X 10-^* 
(0.0061 ± 0.0004) X lO-''* 
(0.58 ± 0.40) X 10-^* 
(0.85 ±0.21) X 10-^* 
(0.06 ± 0.05) X 10-^* 
(0.15 ±0.06) X 10-4* 



where the subscripts indicate which final-state meson contains the spectator quark in the D meson. 
In the diagrammatic approach, amplitudes of the same topology by different subscripts are a priori 
independent of each other. In the SU(3) symmetry limit, however, they are identical. It is also 
instructive to see that in this case, the VP contribution to y in each SU(3) multiplet also vanishes. 
In [l9[, we obtain two possible solutions, called (S,S1) and (A,A1), for the Ty^p, Cv,p, and Ey^p 



amplitudes, depending on which formula is used to extract the invariant amplitudes.^ The ones 
quoted in Table H] are from the (A,A1) solution (upper row) and the (S,S1) solution (lower row). 
We find 

_ f (2.847 ±0.112)% -(2.578 ±0.227)% = (0.269 ±0.253)% (A,A1) 
^^^ ~ j (2.916 ± 0.073)% - (2.764 ± 0.207)% = (0.152 ± 0.220)% (S,S1) ^ ^ 

where the error bars are of the same order as the central values. As far as the central value is 
concerned, yvp is indeed larger than ypp. 



c. vv 

Just as the PP modes, there are 16 VV channels: 4 for CF: K*~p~^,K p ,K lo,K (p; 8 for 
SCS: K*+K*-,K*°K*^,p+p-,p^p^,p^uj,p^(l),ujuj,uj(l} and 4 for DSC: K*+p- ,K*^p^,K*^uj,K*^(j). 
Among them, K (p and K*^cl) are kinematically forbidden, but allowed through the finite width 
effect of K*. Indeed, the decay D^ -^ K (j) has been observed by FOCUS in the Dalitz plot 

^, ZV'+ Ic'- ZV'-^+ 



analysis of D^ -^ K+K-R-tt 

The measurements of i^*^p^,i^ p ,K a; were performed in the early 90 's. During the period 
of 2003-2007, FOCUS had measured K*^ p°4>, K*^K*^ and p°p° through the Dalitz plot analysis 



of various four-body decay modes [25|, |2g]. Some of the VV data are problematic. Naively, it 
is expected from the factorization hypothesis that longitudinal and transverse polarizations are 



comparable in D — )• VV . However, the Mark III measurement 27| has indicated that the branching 
fraction of the D^ — )• K*^p^ decay is already saturated by the transverse polarization state (see, 
e.g., |2S] for a detailed discussion). Since the presently available data do not allow us to have a 
sensible determination of yvvi we will rely on the factorization model to estimate its magnitude. 
Note that the polarized decay amplitudes can be expressed in several different but equivalent 
bases. For example, the helicity amplitudes can be related to the spin amplitudes in the transversity 
basis (^40, ^11 , A±) defined in terms of the linear polarization of the vector mesons, or to the partial- 
wave amplitudes (5, P, D) via: 

1 „ [2 






The solution (S,S1) is obtained by using the equation 



8nni - ^-^ 
to extract the invariant amplitude, while the relation 



87rmi, 

^ pol. 



nD^VP) = ^^\A\\ 



Pc I ;f|2 
87rm? 



is used to get the solution (A,Al), where the polarization vector is taken out of the amplitude so that 
A = {mv/mD)A (e • po)- 

8 



TABLE II: Branching fractions of D^ 
are taken from 



VV calculated in the factorization approach. Data 



151]. Since the W-exchange contributions are neglected in naive factorization, no 

K*^(j) and p^u is made here. 



estimate of the branching fractions of X*''iir ,K ^ sr*o^ „„^ ^o 



Mode 



S wave 
10., 
1 

1.5% 
6.8 X 10-3 

5.8 X 10^3 
0.85 X 10-3 

6.3 X 10-^ 

5.9 X 10-^ 

6.3 X 10-"^ 

3.2 X 10"^ 

4.4 X 10-5 

3.3 X 10-5 



P wave 
6.5 X 10-3 

1.4 X 10-3 
1.2 X 10-3 

5.5 X 10-^ 
6.2 X 10-"^ 

0.91 X 10-^ 
2.5 X 10-5 
6.5 X 10-5 
2.5 X 10-5 
2.8 X 10-5 
2.7 X 10-*^ 
2.0 X 10-6 



D wave 



^total 



Expt. 



K*-p+ 

K*+K* 

p+p- 
p'^pO 

pV 

OJOJ 
UJ(j) 

K*+p- 



1.6 X 10-3 
2.3 X 10-^ 

1.8 X 10-^ 

8.9 X 10-s 

2.3 X 10-'' 

3.4 X 10-5 
1.0 X 10-6 
1.9 X 10-5 
1.0 X 10-6 

4.5 X 10-6 

6.7 X 10-^ 
4.3 X 10-^ 



11. 

1. 

1.6% 
7.3 X 10-3 
6.6 X 10-3 
0.97 X 10-3 
6.6 X 10-^ 
6.8 X 10-^ 

6.6 X 10-^ 
3.5 X 10-^ 

4.7 X 10-5 
3.5 X 10-5 



(6.5 ±2.5)% 

(1.59 ±0.35)% 

(1.1 ±0.5)% 



(1.83 ±0.13) X 10 



-3 



^1 



%/2 



{H+ -H_) = P, 



where we have followed the sign convention of 29[]. The decay rate reads 

Pc 



r{D -^ ^1^2) 



87rm 

Pc 



STrm'jj 



D 

:i^oi'±ivi±p±iviiiP) 



Sttiti'Jj 



(|Sp + |P|2 + |D|2) . 



The factorizable matrix element for the D ^ V1V2 decay is 



^(Wi,y.) ^ (^V2\J^^\o)(^v^\j;,\D) = -ifv,m2 



{el-e*2){mD + mvM?^'iml^,) 



2AP^^ (ml ) ,„ « 

[niD -T mvi ) 
where use of the conventional definition for form factors 



2V^^Hm^ 



V2- 
V2. 



{mo + mvi] 



(20) 



(21) 



(22) 



(/i = 0) and transverse {h = ±) components of Xf^ ^' ^ are given by 



30| has been made. The longitudinal 



X, 



iDVi,V2) 




{m 



D 



-^(Dv^yi) 



ifv2 
2mv^ 



mo + mvi 



mvi 



1 + 



\DVr(2 



Ar'{q^)± 



2pc 



-v^'^Hq') 



(23) 



mo J " " ' mc + mvi 
In the factorization framework, we find that |-ffoP ~ l-f^+P > l-f^-P aiid I'S'P > |-P|^ > l-Dp. 
Therefore, the longitudinal polarization fi defined by 



f -^^ 

11 = -^ 



I A |2 



I U |2 



|Ao|2 + M|||2 + M,|2 



\Ho\^ ± \H+\^ + \H. 



(24) 



is expected to be in the vicinity of 0.5 . Indeed, fi = 0.475 it 0.271 was found in D — )• K*~ p'^ [31|. 
This is not the case in tree-dominated charmful or charmless B — )■ VV decays where the longitudinal 
polarization dominates, i.e., \Hq\'^ > \H^\'^ > |-ff_p and fh = ^ — ^{m'y/m?^). However, for the 
D^ -^ K p^ decay, it was found by Mark III [31] that this mode proceeded through the transverse 
polarization, with only a tiny room for the longitudinal polarization. More precisely, B{D^ — )■ 
^*°P°)trasnverse = (1.6±0.6)%, while the total rate is B{D^ -^ K*° p%ot = (1.59±0.35)%. Mark III 
also measured the partial wave branchin g fr actions: (3.1 it 0.6)%, < 3 x 10^'^ and (2.1 it 0.6)%, 



respectively, for the S-, P- and D- waves |23l. l3ll|. The sum of S- and D-wave branching fractions 



already exceeds the total. Hence, the data associated with this mode are problematic. 

The VV states with different partial waves contribute with different CP parties. We have 
^cp(^^) = 1 for VV in a relative S oi D wave, and —1 in a P wave [14l ]. The parameter y for 
VV modes has the expression 

yvv/ = B{p+p-)t + B{p°py + e(pO^), + ^(/>V), + B{ujuj)e + B{u(f>), + B{K*+K*-), 



+BiK*^K*% -2 cos 6k.+ .- J B{K*-p+),BiK*+p-), - 2cos(5^,o„o VS(i?*V)^^(^*V°)^ 



-2 cos 5k*o^^B{K*Ou;),B{K*Ouj), - 2 cos 5 k*<^ ^^ B{K*y),B{K*y), (25) 

for i = S,D, and an overall minus sign is needed for I = P. 

Using the effective Wilson coefficients ai = 1.22 and 02 = —0.66 19[, the form factors from the 



covariant light-front quark model [34] and the decay constants 33] 



fp = 216 MeV , Jk* = 220 MeV , /^ = 215 MeV , f^ = 187 MeV , (26) 

the predicted branching fractions of D — t- VV decays for various partial waves within the factoriza- 
tion framework are shown in Table [Hi Since the VK-exchange contributions are neglected in naive 
factorization, no estimate of the branching fractions of K*^K , K </>, K*^(j) ^nd p^uj is made here. 
We obtain 

yS-wave ^ Q og% ^ y^-'"'^^^ = -0.03% , 2/^^""*^^ = 0.007% . (27) 

A cancellation between even-parity and odd-parity final states renders small yw in comparison 
with yppyp. 

D. D^ MP 

Recently we have studied the hadronic D meson decays into a pseudoscalar meson P and an 
even-parity meson M, where M represents a scalar meson S, an axial-vector meson A, or a tensor 
meson T [2y]. The data are inferred from detailed Dalitz plot analyses of three-body or four-body 
decays. Normally one applies the narrow width approximation 

T{D -^ MP -^ P1P2P) = r{D -^ MP)B{M -^ P1P2) (28) 

to extract the branching fractions of the D — )• MP decays. There are two complications though: (i) 
Some decays occur near or slightly above the threshold, for example, D^ — )■ /o(980)7r° followed by 
/o — )■ K^K^ . Since the central values of the /o(980) and ao(980) masses are below the threshold 
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TABLE III: Experimental branching fractions for Cabibbo-favored (upper portion) and singly- 
Cabibbo-suppressed (lower portion) decays of D^ — ;■ SP,AP,TP (see (20'1 for details). 



Mode 


Scxp 


Mode 


^cxp 


Mode 


i^cxp 


foK 

17*^ 


(1.6 ±0.5)% 

(8.0+1^) X 10-3 

(8.2 ± 1.4) X 10-3 

(9.2+1^) X 10-3 


i^f(1270)7r+ 
i^-a^(1260) 


(1.14 ±0.32)% 
(7.9 ±1.1)% 


K;-7r+ 


(2.llJ;2) X 10-3 
(S.Otli) X 10-4 


avrO 


(1.0 ±0.3) X 10-4 
(1.8 ±0.3) X 10-4 


7r-a^(1260) 
i^i+(1270)K^ 


(8.98 ± 0.62) X 10-3 
(8.1 ±1.8) X 10-4 


/2vr0 


(3.4 ±0.4) X 10-4 



for decaying into a pair of charged kaons, one cannot apply the narrow width approximation to 
extract B{D^ -^ /o(980)7rO) ^^.^^ jjO _^ fQ{980)TT^ followed by /o -^ K+R-. (ii) Some states, 
e.g., a oi K are very broad in their widths. The use of the narrow width approximation is not 
justified, and it becomes necessary to take into account the finite width effect of broad resonances. 
For example, we find that the branching fraction of D^ — t- cr7r+ extracted from three-body decays 
is enhanced by a factor of 2, whereas B{D^ — )• /2(1270)i^ ) is reduced by a factor of 4 by finite 



width effects 20|]. The current experimental data for the two-body decays of D^ to SP,AP and 
TP are collected in Tables lllli Evidently, one cannot use these data to predict y at this stage. 

The flavor diagram ap pro ach and the factorization calculation have been undertaken to analyze 
these decay processes in [20|. While factorization works well for the CF D^ — )• SP,AP decays, 
predictions are typically about one order of magnitude smaller than experiment for the other decay 
modes, conceivably due to the negligence of weak annihilation contributions arising from FSIs. The 
D —7- TP measurements poise the biggest problem for theory. Predicted branching fractions based 
on factorization are at least two orders of magnitude smaller than data, even when the decays are 
free of weak annihilation contributions. We cannot find possible sources of rate enhancement for 
D^TP. 

There are D — t- MP decays which are not kinematically allowed but may proceed through the 
final widths of even-parity resonances. Examples are D^ — )• Oq (1450)iC^ and D^ — )• K^ (1400)iir^. 
They are needed to ensure the cancellation for y in the SU(3) limit. Beyond the SU(3) symmetry, 
these channels are slightly above the threshold and proceed via finite widths. In this case, SU(3) 
cancellation may be less effective. 



E. Remarks 

Thus far we have concentrated on physical two-body intermediate states. In principle we should 
also consider many body final states. However, we notice that many of 3-body final states arise 
from SP, VP, TP decays, and 4-body states from VV, AP decays. Empirically, we know that 
non-resonant 3-body and 4-body decays are at most 10% of the multi-body decay rate. 

Summing over the currently available data listed in Tables U llIIl we have B{D^ — )• PP) ~ 10%, 
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i3(L>o ^ VP) ~ 28%, BiD° -^ VV) ~ 10%, B{D^ -^ SP) ~ 4.2%, 6(0° -^ AP) ~ 10%, 
B{D^ — )• TP) ~ 0.3%. Hence, the total branching fraction of two-body hadronic decays is 63%. 
This is about 3/4 of the total hadronic rate, recalling that the semileptonic decays account for 16% 
of the total rate 23] ■ This means that, unlike the case of B mesons, the hadronic charm decays 



are dominated by exclusive two-body processes. Since PP and VP final states account for nearly 
half of the hadronic width of D° and ypp+vp = (0.36 ± 0.26)% or (0.24 ± 0.22)%, it is conceivable 
that when all hadronic states are summed over, one could have y ~ (0.5 — 0.7)%. 

III. MASS DIFFERENCE 



Since the short-distance contribution to x is very small [l|, y, [34 1 , it is natural to turn to the 
long-distance effects given by Eq. ([10]): 



Am = V V 



1 {D^\H^\n){n\Hy,\D ) + {D \H^\n){n\H^\D^ 



N 



rriD - E„ 



(29) 



In principle, one can have one-, two-, three-,... particle intermediate states. Unlike the width 
difference, the intermediate states here can be off-shell. The two-body hadronic intermediate state 
contributions to Am were first considered in [ij, [ij] by computing the self-energy diagram depicted 
in Fig. [TJ As stressed in [ij] , the self-energy diagram has a universal imaginary part in the massless 



limit. A dispersive relation between x and y has been derived in [35[ in the heavy quark limit 



Am = V / dE 

2vr J2m^ 



Ar(g) Mqcd \ 

E-TUD \ E J 



Neglecting CP violation in the decay amplitude, a model independent relation ^ 

y ^ 1 - \q/p\ 
X tan (j) 

was obtained in 



(30) 



(31) 



36|-l38|. 
It is convenient to introduce the self-energy correlator 

U{p^) = i [ d'xe'^P-P°>^D'\T[HUx)H^m\D'') • 



(32) 



An insertion of a complete set of intermediate states and use of the representation for the step 
function 



e{t) 



1 



Auit 



-doj 



27ri J uj — ie 
will relate the mass and width differences to the self-energy correlator when p = p£) : 

M^irnl)] ^^_ Im[n(m2,)] 



Am 



AT 



2m,£) 2mD 

The self-energy correlator respects a dispersion relation 

lm[U{s)]ds 



1 r^ 

n(p^) = - / 

TT Jso 



s — p'' + ie 



(33) 



(34) 



(35) 



^ The exact expression of tan (f> was first obtained in 
in the limit of 1 — \q/p\ ^ 1. 



The approximate formula Eq. (|3Tj) is recovered 
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FIG. 1: Two-particle contribution to the neutral charmed meson mass difference, 
where sq = {mi + 777-2)^ with mi and m,2 being the masses in the loop. Hence, 



2mDTD J so s-p^ 
The absorptive part of n(s) amounts to putting the intermediate states in the self-energy diagram 



on shell. Hence, it is proportional to the decay rate of D^ — t- n. The result is 391] 



^ = T^ E ^CKM(n)r?cp(n) cos 5Jb{D^ ^ n)B{D^ ^ n) -^("^i'"^^'^) , (37) 



with 



I(mi,m2,A) = -P / -^ 2-ds , (38) 



so '5 "^D 



where a cutoff A has been introduced to render the loop integral finite and use has been made of 
the formula for two-body decay rates 

B{D' ^n) = -l^\{D^\H^\n)\' . (39) 

A summation over the polarization states of the vector meson is understood for VP final states. 
As shown in [30], the cutoff scale A ~ (2.0 — 2.2) GeV is not far from ?7i£). Note that the result of 
[l3| is recovered in the zero mass limit of intermediate states. 

From Eq. (I37p and Table HI we obtain 

' (0.028 ±0.003)% r (0.064 ±0.009)% for A = 2.2 GeV , 

xpp = < (0.032 ± 0.005)% , xvp = \ (0.073 ± 0.021)% for A = 2.1 GeV , (40) 

, (0.039 ± 0.008)% i (0.088 ± 0.043)% for A = 2.0 GeV . 

Since we have applied the first equation in footnote H] to derive xyp, only the solution (S,S1) is 
relevant for the determination of this parameter. As a result, xppj^vp is of order 10~^ with the 
uncertainty depending on the cutoff scale. As far as the PP and VP modes are concerned, we 
find that xppj^yp is smaller than ypp+vP- This can be seen by comparing Eq. (j37p with Eq. (|12p . 
We see that x is suppressed by a factor of in, while the factor ?7i£)/(mi,m2, A)/pc is maximal for 
D^ —7- TTvr and of order 2.5 . 

Just as the Kl-Ks mass difference receives contributions from the one-particle intermediate 
states such as vr, 77 and ry', ^ it will be interesting to see the single particle effects on charm 



It is known that the contributions to the K^-Ks mass difference from the octet states tt and 778 cancel 
exactly as a consequence of the Gell-Mann-Okubo mass relation, as first noticed in the context of SU(3) 
chiral perturbation theory (ChPT) [40|. A generalization to U(3) ChPT to include the 773 — 770 mixing 



effects was discussed in 



4l| . It turns out that the 77' one-particle intermediate state is one of the main 
contributions to the K^-Ks mass difference, besides the short-distance contribution and the long-distance 
contribution represented by the tt loop [42] ■ 
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mixing. A unique feature of the charm system is that an abundant spectrum of resonances is 
known to exist at energies close to the mass of the charmed meson. Indeed, the sizable magnitude 
of the topological PF-exchange diagram and its large strong phase determined from experiment are 
suggestive of nearby resonance effects. Considering a nearby resonance state R with mass tjir and 
width Tji, its contribution to the mass and width differences is 

{D'>\H^\R){R\H^\D'>) ml-mj 






muF ("^n — 'm,p)2 + m^T'i 



{D^\H^\R){R\H^m Tn , ,.^, 

^^ = ^^ f iml-mlr + mlTl ' ""^ ^^^^ 

XR _ m'jj -mj^ 
VR rnoTR 

where t/r is the CP eigenvalue of the resonance R. One needs to know the weak couplings of the D 
meson to the resonance in order to quantify the resonance contributions to the mixing parameters. 
Some crude estimates had been made in [43(. 

IV. CONCLUSIONS 

Motivated by the possibility of observing new physics effects in D^ — D^ mixing, a lot of efforts 
have been put on experimental determination of the x and y parameters in recent years by the 
BaBar, Belle, and CDF Collaborations. Now that the experimental precision on both quantities 
has reached the level of one per mille, it is timely to scrutinize the SM predictions in order to make 
meaningful inferences. The short-distance contributions to x and y in the SM have been found to 
be several orders of magnitude smaller than the observed values. In contrast, long-distance effects 
from exchanges of multiple hadronic particles play a more important role. 

Since the sum of all two-body hadronic modes that are available currently accounts for about 
63% of the -D" decay branching fraction, it is arguable that these channels dominate and provide 
a good estimate of the mixing parameters. Other multi-body hadronic decays are empirically less 
important, particularly when cancellations among them are taken into account. With more and 
better-determined two-body hadronic decay branching fractions over the past few years, we are in 
a better position to take the exclusive approach to evaluate the long-distance effects. We find that 
the primary contributions to these parameters come from the PP and VP modes. 

To reduce model dependence, we directly take available experimental data and employ the 
diagrammatic approach to evaluate the yet-observed decay branching fractions. From the exchanges 
of PP and VP intermediate states, Eq. (fT7|) for ypp, Eq. (fTU]) for yyp and Eq. pUj) for xpp and 
xvp are the main results of this paper. We obtained that y ~ a few xlO~^ and x ~ 10~^, with 
the latter having a mild dependence on a cutoff scale that is assumed to be around 2 GeV. Here 
we have assumed that the relative strong phase between each pair of Cabibbo-favored and doubly 
Cabibbo-suppressed modes is identically zero, thus maximizing the cancellation. This is partly 
justified by the determination of the relative phase between D^ — )• K^tt^ and K^tt^ decays by 



^ Note that our result of xj^ differs from that in 'l4] by a factor of 2. 
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CLEO and a theoretical estimate of the phase. While inclusive analyses generally render x > y, 
our exclusive calculations indicate that x is smaller than y, in good agreement with the latest direct 
measurements of D^-D^ mixing parameters from the Dalitz plot analysis of D^ — )• Kst^^t^~ and 
KsK~^K~ decays by BaBar. 
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Appendix A: Normalization of states 



In the literature, the D^-D matrix element is usually taken to be 131.135 1 

V 2 A2 2mz)^ ' '"' ^ 2mBV ruD-En + ie ^ ' 

or sometimes ^] 



2/12 ' ' " ' ^ ^h-El + ie 

with different normalizations of the D meson state. It is easily seen that these expressions are 
dimensionally inconsistent. The original formula given by Marshak, Riazuddin and Ryan (MRR) 



m 



M-'-v] =(Z;°|ff^|^°)+y^^°'^-'^^^"'^-'^^ (A3) 

2/12 n mo- En + te 

is dimensionally consistent for one-particle intermediate states. In the MRR convention, the mass 
dimension of the meson is [(/>] = —3/2. In the usual normalization convention with [(/)] = —1, it 
is tempting to put factors of l/2mD on Eq. ()A3p to arrive at Eq. (jAip . However, an appropriate 
normalization factor for the intermediate state \n) {e.g., l/2En for one-particle intermediate states) 
is obviously missing in Eq. (jAip . 

It is instructive to see the derivation of Eq. ()A3p . In the MRR convention, the self-energy 
operator is defined by 

n(p2) = i2mD /'ci^xe*(P-P«)"(DO|r[F^(x)F^(0)]|D°) . (A4) 

It has a canonical mass dimension of 2. An insertion of a complete set of one-particle intermediate 
states 
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and use of the representation ([33]) for the step function lead to 



2m£) 



E/A 



Er, 



6^{qn-pD) 



{D^\H^\n){n\H^\D^) 



+S^{PD -Qn 



Ed- En + ie 

{D'^\H^\n){n\H^\D''y 
Ed - En + ie 



In the rest frame of the D meson, pD = and Ed = inriD- Therefore, 



2mD 



E 



{D'^\H^\n){n\H^\D'^) , {D'^\H^\n){n\H^\D'^) 



■niD - En + ie 



+ 



ruD -En + ie 



(A6) 



(A7) 



Eq. ()A3p then fohows. Since the insertion of (|A5p is dimensionless, this renders Eq. (|A3p dimen- 
sionally consistent. For other normalization conventions and for multiparticle intermediate states, 
an appropriate normalization factor should be included. 



Appendix B: Strong phase of D^ — )■ K^tt relative to D^ -^ K ir^ 

In this appendix we give an estimate of the strong phase S^+n- of D^ — )■ K~^7t~ relative to 
D^ —7- K^TT~^ from the experimental measurement of the former. Prom Eq. (jlSp we learn that 
the strong phase S^+n- vanishes in the SU(3) limit as the double-primed amplitudes T" and E" 
should be the same as unprimed amplitudes T and E, respectively. In the limit of SU(3) symmetry, 
the prediction B{D^ — )• K~^7r~) = (1.12 zt 0.05) x 10~^ is slightly smaller than the experimental 
result of (1.48 ±0.07) x 10~^ [23]. This can be understood as the SU(3) breaking effect in the tree 
amplitude T". In the factorization approach, the relevant tree amplitudes read 

„ Gf 

T = ^( 

\/2 
T" = ^a^ fKiml - mi)F^%mi,) . (Bl) 



:ai/.(m?5-m|)Fo^^(m2) 



^a,fK{ml-ml)F^^{ml 



Taking the form factors for D to vr and K transitions from the recent CLEO-c measurements of 
D meson semileptonic decays to vr and K mesons |4J], we find T" /T = 1.23 . From CF D — t- PP 
decays we obtain (in units of 10~^ GeV) |19l ] 



r = 3.14 ±0.06 



E = (1.53 



+0.07\ „j(122±2)° 
-0.08^ ^ 



(B2) 



Combining the above information, we find that the data of B{D^ — )• i^+vr ) can be better fitted by 
having E" ^ Ee'^^° . This leads to 5K+n~ = arg[(r" ± E")/{T + E)] ^ -Tl^and cos 5K+n~ ~ 0.99, 
consistent with the CLEO measurement of cos Sj^+j^- 



1.03+[]:f^ ± 0.06 
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